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ABSTRACT 



A Hubble Space Telescope (HST) / Advanced Camera for Surveys (ACS) study 
of the structural properties of 85 luminous and ultraluminous (Ltr > 10 114 L Q ) 
infrared galaxies (LIRGs and ULIRGs) in the Great Observatories All-sky LIRG 
Survey (GOALS) sample is presented. Two-dimensional GALFIT analysis has 
been performed on F814W "/-band" images to decompose each galaxy, as appro- 
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priate, into bulge, disk, central PSF and stellar bar components. The fraction of 
bulge-less disk systems is observed to be higher in LIRGs (35%) than in ULIRGs 
(20%), with the disk+bulge systems making up the dominant fraction of both 
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LIRGs (55%) and ULIRGs (45%). Further, bulge+disk systems are the dominant 
late-stage merger galaxy type and are the dominant type for LIRGs and ULIRGs 
at almost every stage of galaxy-galaxy nuclear separation. The mean /-band 
host absolute magnitude of the GOALS galaxies is — 22.64±0.62 mag (1.81q 4 
L}), and the mean bulge absolute magnitude in GOALS galaxies is about 1.1 
magnitude fainter than the mean host magnitude. Almost all ULIRGs have 
bulge magnitudes at the high end (—20.6 to —23.5 mag) of the GOALS bulge 
magnitude range. Mass ratios in the GOALS binary systems are consistent with 
most of the galaxies being the result of major mergers, and an examination of the 
residual-to-host intensity ratios in GOALS binary systems suggests that smaller 
companions suffer more tidal distortion than the larger companions. We find ap- 
proximately twice as many bars in GOALS disk+bulge systems (32.8%) than in 
pure-disk mergers (15.9%) but most of the disk+bulge systems that contain bars 
are disk-dominated with small bulges. The bar-to-host intensity ratio, bar half- 
light radius, and bar ellipticity in GOALS galaxies are similar to those found in 
nearby spiral galaxies. The fraction of stellar bars decreases toward later merger 
stages and smaller nuclear separations, indicating that bars are destroyed as the 
merger advances. In contrast, the fraction of nuclear PSFs increases towards 
later merger stages and is highest in late-stage systems with a single nucleus. 
Thus, light from an AGN or compact nuclear star cluster is more visible at I- 
band as ULIRGs enter their latter stages of evolution. Finally, both GOALS 
elliptical hosts and nearby SDSS ellipticals occupy the same part of the surface 
brightness vs. half-light radius plot (i.e., the "Kormendy Relation") and have 
similar slopes, consistent with the possibility that the GOALS galaxies belong to 
the same parent population as the SDSS ellipticals. 

Subject headings: galaxies: active - galaxies: interactions - galaxies: quasar - 
galaxies: starburst - infrared: galaxies 



1. Introduction 

The Great Observatories All-sky LIRG Survey (GOALS) combines imaging and spec- 
troscopic data from NASA's Spitzer, HST, Chandra and GALEX space-borne observatories 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract No. NAS5-26555. 
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in a comprehensive study of the most luminous infrared-selected galaxies in the local Uni- 
verse (Armus et al. 2009). The sample consists of 181 Luminous Infrared Galaxies (LIRGs, 
L IR = 10 110 - 1O 1L99 L ) and 21 ultraluminous infrared galaxies (ULIRGs, L m > 1O 12 L ) 
derived from the IRAS Revised Bright Galaxy Sample (RBGS, Sanders et al. 2003). The 
sample spans the full range of optical nuclear spectral types (starbursts, LINERs, type-1 and 
type-2 Seyferts; i.e., Veilleux et al. 1995) as well as interaction stages, and serves as a statis- 
tically complete sample of infrared-luminous, local galaxies. As such, the GOALS galaxies 
are excellent analogs for comparisons with infrared and submillimeter-selected galaxies at 
high redshift. 

The Digitized Sky Survey (DSS) images of the RBGS (U)LIRGs at L m > 1O 11O L 
(Sanders et al. 2003) show that more than 90% of them have signs of tidal interaction. 
The galaxy interaction efficiently drives gas to the central regions where vigorous starburst 
activity (Barnes & Hernquist 1992), and possibly the accretion of gas unto supermassive 
nuclear black holes, can commence. Galactic outflows have been found in many (U)LIRGs 
(Armus et al. 1989; 1990, Heckman et al. 1990, Rupke et al. 2005b, Sturm et al. 2011, 
Chung et al. 2011). The galactic scale outflows from stellar winds and supernovae explosions 
play an important role in galaxy evolution. The galactic superwinds clear out nuclear gas 
and dust, enrich the intergalactic medium, and eventually quench star formation activity 
and blackhole growth, and turn gas-rich spiral progenitors into elliptical galaxies or SO-type 
galaxies (Barnes 2002, Veilleux et al. 2005, Sturm et al. 2011) 

Recently, we have undertaken morphological studies of 73 GOALS (U)LIRGs with the 
HST NICMOS H-band images (Haan et al. 2011, hereafter Paper I). We find that a significant 
fraction of the GOALS (U)LIRGs have double (63%) or triple nuclei (6%) which were not 
seen in the B and I band HST data due to obscuration, and the bulge luminosity surface 
density increases significantly along the merger sequence, while the bulge luminosity shows 
a small increase toward late merger stages. This increase in the luminosity surface density 
was found by Haan et al. to be almost entirely due to a decrease of the bulge radius in 
the late stage merging (U) LIRGs, consistent with an inside-out growth of the bulge due to 
the funneling of gas towards the centers. Haan et al. also found that the projected nuclear 
separation is significantly smaller for ULIRGs (median value of 1.2 kpc) than for LIRGs 
(median value of 6.7 kpc), suggesting that the LIRG phase appears at an earlier merger 
stage than the ULIRG phase. The GOALS H-band images are better suited for finding 
multiple nuclei, identifying embedded star clusters and measuring the structural parameters 
in the nuclear region in these dusty GOALS galaxies, but the fields of view are generally too 
small to study the large scale morphologies, bulge-to-disk ratios, tidal features, and stellar 
bars, except in the more distant sources. To address these latter issues, we have undertaken 
morphological studies of the GOALS (U)LIRGs with the HST ACS I-band images. Optical 
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images remain excellent tools for measuring bar fractions and bar length, as long as the 
rest-frame wavelengths are long wards of the Balmer break (e.g., Sheth et al. 2003, 2008; 
Menendez-Delmestre et al. 2007, 2012). The present paper aims to study optical structural 
properties of GOALS (U)LIRGs with L m > 10 114 L using GALFIT (Peng et al. 2002, 
2010). 

The ACS on HST provides superior spatial resolution, making it possible to image 
detailed structural features such as unresolved nuclei, faint tails, bridges, rings, and shells. 
In addition, the large field of view (202" x 202") of the Wide Field Channel (WFC) on 
ACS makes it possible to image (U)LIRGs consisting of widely separated pairs in a single 
orbit. The structural decomposition of GOALS (U)LIRGs allow an assessment of the relative 
fractions of (U)LIRGs with disks, disks+bulges, and elliptical-like profiles, the structural 
properties of the progenitors (early stage mergers) and evolutionary byproducts (late-stage 
mergers) of LIRGs and ULIRGs, and their fundamental plane properties. Further, the 
present sample spans all the merger evolutionary stages (e.g., Sanders et al. 2003), making 
it possible to assess these properties as a function of merger phase. Note, however, that our 
sample was selected, in part, by luminosity cut (L IR > 10 1L4 L ) and thus may be biased 
toward later stage mergers. The results presented here are representative of all (U) LIRGs 
with Lir > 10 114 L in the local universe. 

The paper is divided into five Sections. In §2, the sample selection and observations are 
summarized. In §3, the data reduction and analysis are described. The results and discussion 
are presented in §4 & §5, and a summary of the paper is presented in §6. Throughout this 
paper, the cosmology H = 70 km s _1 Mpc" 1 , VLm = 0.3, and = 0.7 are adopted (see 
also Armus et al. 2009). 



2. Sample Selection and Observations 

Sample galaxies were selected from a complete sample of 87 (U)LIRGs with Lj^ > 10 114 
L Q in the IRAS Revised Bright Galaxy Sample (RBGS, i.e., > 5.24 Jy and Galactic 
Latitude |6| > 5°, Sanders et al. 2003). These (U)LIRGs have been imaged with HST with 
the ACS/WFC using the F435W (B) and F814W (I) broad-band filters (GO program 10592, 
PI: A. Evans: see Evans et al. 2013). One galaxy was observed per orbit in the ACCUM 
mode, with total exposure times of ~ 21 minutes and ~ 12 minutes in F435W and F814W 
filters, respectively. The F435W and F814W observations were done using the three and 
two point line dither patterns, respectively. Further details of the observations can be found 
in Evans et al. (2013). The basic properties of GOALS galaxies in our sample are listed in 
Table 1. 
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3. Data Reduction and Analysis 

The data reduction was initiated with the flat-fielded science file retrieved from the 
Multimission Archive at Space Telescope (MAST). The four quadrants of the flat-fielded 
images showed a varying range of bias level offsets (up to 20%); these quadrant offsets were 
measured and corrected, and the sky background was set to zero. Next, the cosmic-rays in the 
images were removed with the lacosmic routine (van Dokkum 2001), which uses a Laplacian 
edge detection algorithm. The algorithm works better on non-drizzled images since drizzling 
smooths the sharp edges. Individual images with the cosmic-ray removed were combined 
using Multidrizzle package (Koekemoer et al. 2002); this routine also removes additional 
cosmic-rays and corrects the geometric distortion of the images. The remaining cosmic rays 
were further processed with jcrrej2 routine in IRAF (Rhoads 2000). Finally, the combined 
images were rotated and WCS header informations were corrected by making use of 2MASS 
catalog coordinates of bright stars in each image. 

Model fitting of the images was done using the 2-dimensional galaxy fitting package 
GALFIT. The procedure was used to fit the central point source in each object and deter- 
mine the structural parameters of the underlying host galaxies. The analysis of each object 
followed a number of well defined steps. First, a mask was constructed to exclude bright 
stars or small foreground/background galaxies within the field of view. Next, because we 
are particularily interested in disk and bulge properties in underlying host galaxies, a fit was 
made to the surface brightness profile of each object for the following cases of the disk and 
bulge combinations: 

(%) PSF + exponential disk profile (n = 1) alone, 
(%%) PSF + de Vaucouleurs profile (n = 4) alone, and 
(Hi) PSF + disk (n = 1) and bulge (n = 4), 

where PSF is the point-spread-function and n is the power-law index of the Sersic profile 
(Sersic 1963), i.e., 

S(r) = S e e^ r /^ 1/n - 1 ]. (1) 

In Equation (1), r e and S e are the effective radius and surface brightness, respectively. In 
cases where a bar can visually be identified in a disk and disk+bulge galaxies (i.e., 25.9% of 
the sample), a bar component was included in the 2-D fit. Next, chi-square values of the fit, 
fitted components, and residual image for the above 3 cases were examined to identify the 
best fit model among these 3 cases. For the galaxies with multiple nuclei, the fitting was 
done simultaneously. 
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Early in the fitting process, it was determined that the fits to the F435W images did 
not converge; the dusty nature of these (U)LIRGs renders their appearance at 0.4/xm too 
patchy and discontinuous for a GALFIT analysis. These fits were thus abandoned - only 
fits to the F814W data are included and discussed here. Note also that nuclei of two of 
the 87 (U)LIRGs observed with HST, i.e., Mrk 231 and IRAS 05223+1908, are saturated 
at F814W, and thus reliable fits for these (U)LIRGs could not be achieved. Fig. 1 contains 
an example image of the galaxy, the model version of the galaxy, the residual (i.e., galaxy - 
model fit), and the individual model components (PSF, bulge, disk, and bar). Table 2 lists 
structural parameters derived from the best fit model for each of the remaining 85 (U)LIRGs 
(total number of nuclei is 137: 38 singles + 42 doubles + 5 triples) in the F814W filter. 

4. Results 

As stated earlier, the sample contains wide range of interaction classes from advanced 
mergers to widely separated systems with nuclear separation up to 65 kpc. In this sense, the 
GOALS sample represents an evolutionary sequence for interacting galaxies. The current 
GOALS sample consists of 64 LIRGs (75.3%) and 21 ULIRGs (24.7%) and will be discussed 
in terms of these luminosity classes where appropriate. For the objects with multiple nuclei, 
the infrared luminosity for each nucleus will be apportioned according to its Wide-field 
Infrared Survey Explorer (WISE) 22/xm flux or Two Micron All Sky Survey (2MASS) Ri- 
band magnitude. Throughout the paper, the term "single- nucleus merger (U) LIRGs" is used 
to refer to (U) LIRGs in which the nuclei of their progenitors have coalesced and thus a single 
nucleus is observed. 



4.1. Galaxy Morphological Types 

The fractions of morphological types derived from the GALFIT analysis for the GOALS 
host galaxies are listed in Table 3 and plotted in Fig. 2. The fraction of disk, disk+bulge, 
and elliptical in GOALS galaxies are 32.1%, 48.9%, and 19.0% respectively. If we subdivide 
GOALS sample into LIRGs and ULIRGs, the distributions of morphological types are rather 
different. Both classes have a similar fraction of disk+bulge galaxies, but a higher fraction 
of disks are found in GOALS LIRGs, whereas a higher fraction of ellipticals are found in 
GOALS ULIRGs. The mean (range) projected nuclear separations of the GOALS LIRGs 
and GOALS ULIRGs are 12.2 kpc (0 - 65.1 kpc) and 4.0 kpc (0 - 39.2 kpc) respectively. This 
suggests that more luminous, advanced merger systems tend to have a relatively smaller disk 
fraction and a higher elliptical fraction. Note, however, that the disk+bulge galaxy types 
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make up the higher fraction of both LIRGs (55%) and ULIRGs (45%). 

The distribution of morphological type as a function of interaction class (see Haan et 
al. 2011; Evans et al. 2013 for more details) is plotted in Fig. 3. The following is a brief 
description of interaction class (IC) based solely on a visual inspection of each systemj^] 

- single undisturbed galaxy, shows no signs of tidal interaction 

1 - separate galaxies, disks symmetric (intact), no tails 

2 - progenitor galaxies distinguishable, disks asymmetric or amorphous, tidal tails 

3 - two nuclei in a common envelope 

4 - double nuclei + tidal tail 

5 - single or obscured nucleus, long prominent tails 

6 - single or obscured nucleus, disturbed central morphology, short faint tails or tails 
absent, shells. 

At the widely separated stage (IC=1), elliptical hosts occupy the smallest fraction of 
morphological types (11%) whereas disk+bulge systems occupy the largest fraction (61%). 
The fraction of disk galaxies gradually increases from IC=1 to IC=3 and decreases from IC=3 
to IC=5. In contrast, the disk+bulge fraction decreases from IC=1 to IC=3 and increases 
from IC=3 to IC=6, and, with the exception of IC = 3, they make up the highest fraction 
of (U)LIRGs at every interaction stage. In general, there is no clear trend in morpholog- 
ical types as a function interaction class, but we find more elliptical fraction in advanced 
single-nucleus merger (IC=5) than in widely separated galaxies (IC=1). The distribution of 
morphological types in isolated systems (IC=0) is similar to that in advanced mergers (IC=5 
& IC=6) and all but one galaxy in IC=0 are either disk+bulge or elliptical galaxies. This 
may be due to the fact that many IC=0 systems are late-stage mergers in which the tidal 
features have faded. 



4.2. Host Properties 

Fig. 4 shows distributions of host (total - PSF component, left panel) and bulge (right 
panel) absolute magnitudes for GOALS galaxies (red), GOALS LIRGs (green), and GOALS 
ULIRGs (blue) respectively. The vertical dashed line represents M\ = —22.0 mag (Blanton 
et al. 2003), i.e., the J-band absolute magnitude of an L* galaxy in a Schechter luminosity 



2 More detailed modeling of the interaction stages which takes into account kinematics is underway 
(Privon, Ph. D. Thesis, Privon et al. 2013). This modeling will allow us to address obvious degeneracies in 
our present classification scheme. 
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function of the local galaxies. The host absolute magnitudes (for double and triple galaxies, 
all components were summed since these interacting components will eventually be com- 
pletely merged to form a single nucleus system) for the GOALS galaxies range from —20.98 
to —23.94 magnitudes, with a mean (median) of — 22.64±0.62 (—22.75) mag. There is es- 
sentially no difference in the mean absolute magnitudes between the GOALS LIRGs and 
GOALS ULIRGs (-22.61±0.68 mag for LIRGs and -22.72±0.41 mag for ULIRGs). The 
average magnitude of the GOALS host galaxies corresponds to 1.8lJ;4 Lj- 

The bulge absolute magnitudes of the disk+bulge and elliptical hosts range from —17.54 
to —23.70 mag with a mean (median) of — 21.56±1.03 (—21.77) mag. The mean bulge lumi- 
nosity is about 1.1 mag fainter than the mean host galaxy luminosity. Although the distribu- 
tions are different, we do not find a significant difference of mean bulge absolute magnitude 
between the GOALS LIRGs and GOALS ULIRGs (-21.38±1.06 and -22.07±0.78 for LIRGs 
and ULIRGs respectively). 

The distribution of bulge-to-disk ratio, B/D, is plotted in Fig. 5. The B/D ratios (i.e. 
ratio of Sersic index of n=4 bulge component to n=l disk component from GALFIT) range 
from 0.01 to 7.9 with a mean (median) of 1.23±1.66 (0.67). Most of the B/D in GOALS 
galaxies are less than 2, and more than 60% are less than 1, thus the majority of the GOALS 
galaxies are disk-dominated system. The median B/D value of 0.67 in GOALS galaxies is 
comparable to that of the local Sa-Sab galaxies (Oohama et al. 2009). 

In binary systems, a mass ratio between two interacting galaxies can be estimated from 
their luminosity ratio by assuming that they have the same M/L ratio. Fig. 6 shows the 
distribution of mass ratio in the binary systems in interacting galaxies (Fig. 6a) and the 
mass ratio as a function of M(host) (Fig. 6b). The result suggests that most of the GOALS 
binaries are major mergers (i.e., mass ratio of fTi secon( iary/ ^primary ^1/4) and only 7% are 
minor mergers. The mean value of the mass ratio is 0.55±0.23 in all GOALS binaries. The 
mean mass ratio of the GOALS ULIRGs (0.60±0.23) is similar to that of the GOALS LIRGs 
(0.55±0.23). In general, the absolute magnitude of the equal mass interacting galaxies is 
brighter than that of the unequal mass interacting galaxies. However, this correlation is very 
weak (r=0.34) as shown in M(host) and mass ratio plot (Fig. 6b). 



4.3. Stellar Bars 

Bars are ubiquitous in disk galaxies and play an important role in the dynamical and 
secular evolution of galaxies. A series of GALFIT measurements with a bar component 
included was run on the /-band images of the GOALS sample. GALFIT successfully de- 
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composed bars for the isolated and mildly disturbed systems, however, such a decomposition 
was unsuccessful for very disturbed systems. For these disturbed systems, we have relied on 
visual inspections to identify the bars. The result of bar decomposition is listed in Table 2: 
if the bars were identified bar by visual inspection, it is flagged as 'Y' in column 9 and if 
not, it is flagged as 'N'. If the bar components were successfully fitted with GALFIT, it is 
flagged as 'B' in column 2. 

4-3.1. Bar fraction 

The bar fraction in nearby, bright spiral galaxies is ~65% (de Vaucouleurs 1963) with a 
third classified as strong (SB) bars and another third classified as intermediate (SAB) bars. 
A variety of studies have confirmed this fraction (e.g., Moles et al. 1995; Mulchaey & Regan 
1997; Hunt & Malkan 1999; Eskridge et al. 2000; Knapen et al. 2000; Laurikainen et al. 2004; 
Sheth et al. 2008). Studies in the infrared, where the stellar bars are more easily identified, 
have also confirmed this bar fraction and have shown that many of the "intermediate" (SAB) 
bars in the optical are re-classified as "strong" (SB) bars in the infrared (e.g. Eskridge et 
al. 2000; Menendez-Delmestre et al. 2007; Sheth et al. 2012) but the overall bar fraction in 
bright, nearby disk galaxies is unchanged. At the same time, some studies which have quoted 
lower bar fractions in spirals (e.g. 48-55%, Aguerri et al. 2009 or Barazza et al. 2008) can 
be attributed to different, and less accurate bar identification techniques or different sample 
selection criteria. In comparison to these largely non-interacting, normal disk galaxies, the 
overall bar fraction in the GOALS sample is much lower (25.9% - 29 out of 112 galaxies). 

The bar fraction as a function of interaction class and morphological type are presented 
in Fig. 7 (fraction on top of each histogram represents number of galaxies with bar to the 
total number of galaxies in the bin) and tabulated in Table 4. The bar fraction gradually 
decreases from widely separated systems to more advanced merger systems. When galaxies 
are relatively well separated (IC=1), about 56% of the galaxies have bars. This fraction is 
similar to that in nearby spiral galaxies. However, if the galaxy orientation is considered 
(see below), this bar fraction represents a lower limit. In most cases, the studies of bar 
fraction have been restricted to galaxies having inclination angles less than i < 60 ~ 70°. If 
the GOALS galaxies are randomly orientated and a reliable bar decomposition can only be 
applied to galaxies with orientation angle < 70°, then the bar fraction of 56% found in well 
separated galaxies (IC=1, mean nuclear separation=~ 40 kpc) would effectively be boosted 
to 72% (9/7x56%). This fraction is about 10% to 20% more than that found in optical 
images of spirals; if the simplistic extrapolation is correct, it is suggestive of tidally-induced 
bar formation (Toomre & Toomre 1972; Noguchi 1987; Gerin et al. 1990; Steinmetz & 
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Navarro 2002) already commencing at nuclear separation of ~ 40 kpc. 

The mean bar fraction in completely merged systems (IC=5 and IC=6) is much smaller 
than that of the whole sample suggesting possible destruction of bars at the terminal stage 
of the interaction. The bar fraction in IC=0 is 28.6% which is approximately 5 times larger 
than that of the single nucleus mergers (IC=5 and IC=6) but about a half of that found 
in nearby spirals. In Fig. 7b, the bar fraction in GOALS disk+bulge systems (32.8%) is 
about twice that observed in disk galaxies (15.9%). If we calculate B/D ratios in bulge+disk 
systems that contain bar, most of them are disk-dominated systems (B/D <1) and their 
median value of B/D is only 0.18. This value is about 1/3 of B/D value found in total 
disk+bulge systems, suggesting bar in disk+bulge systems are preferably found in GOALS 
galaxies with a large disk and a small bulge. 



4-3.2. Strength, size, and ellipticity of bars 

We were able to fit bar components for 19 GOALS nuclei with GALFIT. The distribution 
of bar-to-host intensity ratio Ibar/^-host is shown in Fig. 8a. The Ibar/^host ranges from 0.4% 
to 30.2% with mean value of 12.8+7.8%. Similar values are found in nearby field galaxies 
by Reese et al. (2007, 17.7+11.3%) and Weinzirl et al. (2009, 16.6+10.2%). Fig. 8b shows 
the distribution of bar half-light radius r e (bar). The bar half-light radius in the sample 
ranges from 0.3 kpc to 6.1 kpc with mean of 2.7+1.6 kpc, consistent with other studies of 
bar sizes (e.g., Menendez-Delmestre et al. 2007; Weinzirl et al. 2009). The similar values 
of Ibar/^-host and r e (bar) in GOALS galaxies with those found in nearby spirals suggest that 
bars in galaxies have characteristic ranges of har/^-host and r e (bar) regardless of the source of 
perturbation. 

Fig. 9 shows the distribution of bar ellipticity e (Fig. 9a), r e vs. e (Fig. 9b), and 
B/D vs. e (Fig. 9c) respectively. The bar ellipticities in our sample range from 0.36 to 
0.88 with an average value of 0.68+0.14. The GOALS galaxies value is consistent with what 
prior studies have found in nearby galaxies surveys: e=0.67±0.14 (Weinzirl et al 2009), 
e=0.50±0.13 (Menendez-Delmestre et al. 2007), e=0.50±0.10 (Marinova & Jogee 2007), 
and e=0.49±0.14 (Reese et al. 2007). It is interesting to note that bar ellipticities derived 
from GALFIT (this paper & Weinzirl et al. 2009) are about 30% larger than those derived 
from other methods (nonparametric decomposition method, Reese et al. 2007; and ellipse 
fit method, Menendez-Delmestre et al. 2007; Marinova & Jogee 2007). This may reflect the 
inability of non-decomposition methods to adequately remove the bulge contribution in the 
measurement of bar ellipticity. 
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In Fig. 9b, a statistically significant correlation (r=0.61) is found between the bar 
ellipticity and bar half-light radius. The small bars tend to be round and the large bars tend 
to be elongated. The rarity of ellipticities less than 0.4 suggests that either it is difficult 
to form (or detect) ovals in merging galaxies - most bars in merging systems are relatively 
strong. In Fig. 9c, a weak correlation (r=0.36) is observed between the bar ellipticity 
and B/D, i.e. bars become round (small ellipticity) in galaxies having large bulges. This is 
consistent with Das et al. (2003) but we note that this is likely an observational bias because 
larger bulges lead to apparently fatter bars (e.g., see for example Figure 2 in Sheth et al. 
2000). 

In order to find any distinctions that may exist between the barred and unbarred sys- 
tems, the mean values of M(host), M(bulge), and B/D have been estimated for both popula- 
tions (Table 5). The mean M(host) in barred systems is similar to that in unbarred systems 
(Am = 0.23 mag). On the other hand, the mean M(bulge) in barred systems is about one 
magnitude fainter (Am = 0.78 mag) than that in unbarred systems. This may indicate the 
stabilizing influence of large bulges in merging galaxies against bar formation (e.g., Mihos 
& Hernquist 1996). However note that the frequency of bars in early type galaxies (more 
bulge-dominated systems) is higher at higher redshifts (z~l) than late type galaxies (Sheth 
et al. 2008) but in the local Universe, the fraction of bars in bright, disk galaxies is roughly 
the same across Sa-Sc galaxies (Eskridge et al. 2002; Menendez-Delmestre et al. 2007 and 
references therein; Sheth et al. 2012). The role of bulges in formation or inhibition of bars 
during interactions and mergers is worth exploring in more detail with larger samples. 



4.4. Nuclear Point Source 

As galaxies reach more advanced merger stages, more gas will be driven to the central 
region and thereby provide a more favorable environment for compact nuclear starburst and 
AGN activity. In this section, the properties of the point source (PSF) component in the 
GOALS galaxies are examined. The PSF magnitudes derived from GALFIT are listed in 
column 8 of Table 2. The FWHM of the PSF component is about 0.1 arcsec (~2 pixels) 
and mean redshift of the GOALS galaxies is z=0.037. Thus, the maximum physical size of 
the corresponding region of the PSF emission is ~ 70 pc, i.e., much smaller than that of the 
typical compact, nuclear starburst region in (U)LIRGs (a few hundred pc). 
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4-4.1. Detection rate of PSF component 

A PSF component has been detected in 22 out of 137 GOALS nuclei (16.1%). The 
fractions of galaxies with a PSF component in GOALS galaxies as a function of interaction 
class, morphological type, and infrared luminosity (LIRG vs. ULIRG) are plotted in Fig. 
10 and tabulated in Table 6. The fraction on top of each histogram represents number of 
galaxies with PSF detection to the total number of galaxies in the bin. The PSF fraction in 
interacting galaxies (IC=1 to IC=4) ranges from 5.9% to 16.7%. When they become more 
advanced single nucleus mergers, the PSF fraction increases rapidly (33.3% in IC=5 and 
50.0% in IC=6). If we compare PSF fractions in interacting galaxies (IC=1 to IC=4) and 
more advanced single nucleus mergers (IC=5 and IC=6), we can have more robust statistics 
in each group with the increased number of data. The PSF fraction in single nucleus mergers 
(9/25 = 36.0%) is about 4 times larger than that in interacting galaxies (10/103 = 9.7%). 
The PSF fraction in isolated systems (IC=0, 33.3%) is similar to that in single-nucleus 
mergers (note again that the isolated systems might be post-mergers, but the total number 
of galaxies in this bin is not sufficient to confirm; only 3 out of 9). 

The PSF fractions are 2.3% (1/44), 22.1% (15/67), and 25.9% (6/26) for disk, disk+bulge, 
and elliptical galaxies, respectively. Compared to the disk galaxies, the PSF detection frac- 
tions in disk+bulge and elliptical systems are about an order of magnitude larger. If the 
PSF component is a compact nuclear starburst, the above result suggests that elliptical-like 
GOALS galaxies have more compact starbursts or unresolved star clusters in the nucleus 
than in disk-like GOALS galaxies. 

The PSF fraction in ULIRGs (25.0%) is about 1.8 times larger than that in LIRGs 
(13.8%). In general, the molecular hydrogen mass and star formation rate are higher in 
ULIRGs than in LIRGs. Large and more concentrated molecular gas in ULIRGs than in 
LIRGs implies that ULIRGs have more dust extinction than in LIRGs. This suggests that 
outflows associated with nuclear activity in ULIRGs is strong enough to initiate clearing out 
the nuclear gas. 

44.2. PSF strength 

Distributions of PSF to host intensity ratios in GOALS galaxies are shown in Fig. 
11. The IpsF/Ihost in GOALS galaxies estimated from I-band images (left panel) ranges 
from 0.002 to 0.056 with mean (median) of 0.015±0.012 (0.014). The I PSf /I w in GOALS 
elliptical galaxies (cyan line on top-left panel) tend to be larger than the mean value. We do 
not see any difference in the IpsF/hiost between GOALS LIRGs (median=0.015) and GOALS 
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ULIRGs (median=0.013). The right panel in Fig. 11 shows distributions of IpsF/^-host 
estimated from H-band images (Haan et al. 2011) for the GOALS galaxies classified as 
ellipticals (14 objects listed in Table 8). We might expect that the IpsF/^host would be 
larger as measured in the H-band compared to that measured in the I-band simply due to 
the reduced extinction in the near-infrared (extinction at H is 1/3 less than at I). Among 
the GOALS LIRGs we see little or no difference between the I-band and H-band data, but 
the GOALS ULIRGs appear to have a larger ratio as measured in the H-band compared to 
the GOALS LIRGs seen in both filters. The median Ipsr/^host among GOALS ULIRGs is 
0.033 while only 0.008 among GOALS LIRGs. 

Fig. 12a shows mean IpsF/hiost as a function of morphological type. The IpsF^host 
of GOALS elliptical hosts is about twice as large as that of disk+bulge hosts. One of the 
main differences between disk+bulge and elliptical systems is the bulge fraction, and if the 
ellipticals have stronger PSF emission than disk+bulge systems, the PSF emission might be 
correlated with the bulge mass. Fig. 12b shows absolute bulge magnitude versus absolute 
PSF magnitude plot. We find a marginally good correlation (r=0.54) between PSF and 
bulge magnitudes. Correlations have been measured between the black hole mass and bulge 
magnitude and the nuclear star cluster magnitude and bulge magnitude (Ferrarese et al. 
2006); the PSF-bulge magnitude correlation observed for our sample are consistent with 
both. 



4.5. Tidal features 

Tidal features such as tails, bridges, lopsided disks, distorted outer isophotes, and off- 
center nuclei are visible in almost all of the GOALS galaxies. The strength of the tidal fea- 
tures has been estimated from the residual maps in Fig. 1 and the intensity ratio of residual 
to host galaxies is listed in Column 12 of Table 2. Distributions of residual to host intensity 
ratios in GOALS galaxies (red), GOALS LIRGs (green), and GOALS ULIRGs (blue) are 
shown in Fig. 13. The l res /^host in GOALS galaxies ranges from 0.18 to 0.71 with mean 
(median) of 0.38+0.11 (0.36). We do not find any difference in l res /l host between GOALS 
LIRGs (mean (median)=0.37±0.11 (0.35)) and GOALS ULIRGs (mean (median)=0.40±0.12 
(0.38)). These numbers are, of course, consistent with GOALS galaxies being greatly dis- 
turbed systems. 

Fig. 14 shows l res /lhost of each component in GOALS binary systems as a function of 
host absolute magnitude. In this Figure, binary pairs are connected with a solid line. The 
circles with the red bar and the green filled circles represent galaxies with a stellar bar and 
galaxies with a PSF, respectively. In a binary system, if the l res /^host in a small companion 
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(faint component) is larger than that in large one (bright component), it is plotted in Fig. 14a 
and the opposite case is plotted in Fig. 14b. In GOALS binary systems, the fraction of the 
systems where the small companion has a larger \ res /^host (suffering more disturbance) than 
the large one is 60%, and it is 40% for the opposite case. About 30% more bars and about 
twice as many PSFs are found in the binary systems where the small companion suffers more 
disturbance than the large one (galaxies plotted in Fig. 14a). The mean l res /\h ost values of 
the binary systems in Fig. 14a (I res /I ftost =0.41±0.18) and in Fig. 14b (I res /I ftost =0.36±0.15) 
do not show any significant difference. 

5. Discussion 

5.1. Trend With Nuclear Separation 

In this section, we will investigate trend of morphological parameters as a function of 
nuclear separation. For this analysis, we will use single and double nucleus systems (isolated 
systems, IC=0, are not included). The nuclear separation is divided into 4 bins in order to 
have a comparable number of data points (0 kpc, 0-10 kpc, 10-30 kpc, 30-70 kpc). 

We have plotted 9 basic morphological parameters as a function nuclear separation in 
Fig. 15 (fraction of morphological types (a), host magnitude (b), huige/^-host (c), bar fraction 
(d), l bar /l host (e), bar size (f), PSF fraction (g), l PS F/host (h), and l res /I host (i)). 

Simulations suggest that strong starburst activity commences in interacting galaxies 
after their first pericentric encounter (which is around 50 kpc, see details in Gerritsen 1997; 
Mihos & Hernquist 1996). At this distance, the dominant morphological type in GOALS 
galaxies is the disk+bulge system (NS=30-70 kpc in Fig. 15a, disk+bulge systems = ~70%, 
and disk galaxies = ~25%). If the morphological types of widely separated GOALS galaxies 
are representative of (U)LIRG progenitors, then the majority of progenitors of the GOALS 
galaxies are disk+bulge systems. Elliptical hosts represent about 20% of the galaxies when 
the nuclear separation of the GOALS galaxies becomes less than 10 kpc. The absolute 
magnitude of the host galaxies as a function of nuclear separation (Fig. 15b) shows no 
trend. The M(host) in galaxies with small nuclear separation could be smaller than that in 
galaxies with large nuclear separation due to increased extinction. However, we do not find 
any difference of M(host) between galaxies with NS > 10 kpc (M(host)=— 22.51 ±0.59) and 
galaxies with NS < 10 kpc (M(host)=-23.07 ± 0.50). 

Although there exist large error bars in \uigel\ost vs. NS plot (Fig. 15c), we find 
increasing trend in the median value of huige/Ihost with decreasing nuclear separation. The 
median values of Ibuige/^host in each NS bin are 0.18, 0.24, 0.53, and 0.60 for galaxies with 
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NS=30-70 kpc, NS=10-30 kpc, NS=0-10 kpc, and NS=0 kpc, respectively. When galaxies 
are relatively well separated, the huige/^-host increased slightly (~33% from NS=50 kpc to 
NS=20 kpc). However, when they are close together, a significant increase of huige/Ihost 
occurs (250% from NS=20 kpc to NS=0 kpc). If we assume galaxies at NS=20 kpc become 
completely merged (NS=0 kpc) with radial velocity of 200 km s -1 , the huige/^-host doubles in 
7.8 xlO 7 years. 

The bar detection fraction in GOALS galaxies (Fig. 15d) shows a decreasing trend as 
a function of nuclear separation. When galaxies are relatively well separated (NS=30 — 70 
kpc), about 57% of the galaxies have a bar. This fraction is similar to what have been found 
in other samples of nearby spiral galaxies. If we consider orientation effects as discussed in 
§4.3.1, the bar fraction of 57% in well separated galaxies will become 73% (9/7x57%), which 
is about 10% larger than that in normal galaxies. Fig. 15d is possible evidence that the bar 
destruction rate increases with decreasing nuclear separation: ~ 20% of bars are destroyed 
in the first 50 kpc to 20 kpc interval (30 kpc), ~ 25% in the following 15 kpc interval, and 
~ 15% in the final 5 kpc interval. On average, it will take about 1.2xl0 8 yrs to destroy half 
of the bars from NS = 50 kpc NS=0 kpc if we adopt mean interaction speed of 200 km s -1 . 
This value is about 10% of the bar decaying time in Sb and Sc galaxies as estimated from 
numerical simulations (Bournaud et al. 2005). 

If the bars are decaying slowly, we should observe a smooth trend of decreasing har/^-host 
ratio as a function of nuclear separation among the GOALS galaxies that show bars. On the 
other hand, if the bars are decaying abruptly, no dependence of Ibar/^host with decreasing 
nuclear separation will be observed. As seen in Fig. 15e, we do not see any dependence of 
har/^-host on the nuclear separation, suggesting bar destruction could occur rather abruptly. 
Similarly, the fact that we see no correlation of the bar half light radius with nuclear sepa- 
ration argues against a slow decay of existing bars (Fig. 15f). 

The fraction of GOALS galaxies with PSF component as a function of nuclear separa- 
tion is plotted in Fig. 15g and tabulated in Table 6. The PSF detection fraction for widely 
separated galaxies (NS = 30 — 70 kpc) is about 17% whereas it is about 8% for closely 
interacting (NS=10-30 kpc) and overlapping disk galaxies (NS=0-10 kpc). The PSF detec- 
tion fraction is expected to increase with decreasing nuclear separation since the gas that 
fuels vigorous starburst or related AGN activity will be more centrally-concentrated in the 
nuclear region. An additional contributing factor that causes the opposite effect is obscu- 
ration. So, compared to the widely separated galaxies, the small PSF detection fraction in 
closely interacting and overlapping disk galaxies might be caused by increased extinction. 
This can become more apparent if we compare PSF detection fractions of the same objects 
observed in I and H bands. For the GOALS elliptical galaxies, we find a PSF detection 
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fraction of 26.3% with I-band images and 73.7% with H-band images of Haan et al. (2011). 
If a — 1.5 is assumed for the extinction law of Ax oc X~ a , extinction at /-band is about 
3 times larger than that in if-band. We find a sharp rise of the PSF detection fraction 
from overlapping disk galaxies to single- nucleus merged galaxies (8% to 38%, i.e., a 5 times 
increase). The mean IpsF/^-host increases from widely-separated galaxies (NS=30-70 kpc) 
to overlapping disk galaxies (NS=0-10 kpc), but decreases at NS=0 kpc probably because 
of strong dust obscuration in the nuclear region (Fig. 15h). Note, however, that the large 
error bars also allow for the possibility that there is no change in IpsF/^-host as a function of 
nuclear separation. 

A low overall fraction of light in tidal features is expected as a function of decreasing 
nuclear separation. However, we do not see such trend in l res /Ih st vs. nuclear separation 
plot (Fig. 15i), suggesting a significantly long time is required to settling down the tidal 
features even when they are completely merged. A similar result was also observed in 1 Jy 
ULIRGs, where 78.5% (51 out of 65) of completely merged ULIRGs (Veilleux et al. 2002) 
show prominent tidal features. Simulations (Hibbard & Mihos 1995) and deep optical & 
radio observations (Hibbard & van Gorkom 1996) of interacting galaxies suggest that the 
bases of the tails fall back quickly to galaxies, while the more distant regions fall back more 
slowly. This delayed return of tidally ejected material may evolve on very long time scales, it 
may leave observable signatures for many Gyrs (Hibbard & van Gorkom 1996) as evidenced 
in tidal features of (U)LIRGs. The infalling tidal tails could be the source of disk growth 
over long period of time (Barnes 2002). 



5.2. Kormendy Relation 

Fig. 16 shows surface brightness vs. half-light radius for the elliptical hosts of the 
GOALS galaxies (red circles) and SDSS ellipticals (green dots; Bernardi et al. 2003 using 
/ — 2=1.00). The least-square fit of the GOALS ellipticals (red line) shows similar slope as 
that of the z-band data of Bernardi et al. (green line). The mean size of GOALS galaxies is 
ri/ 2 =5.3±3.7 kpc, which is comparable to that of the SDSS ellipticals (ri/ 2 =5.2±3.0 kpc). 
These two results may be an indication that the GOALS galaxies may belong to the same 
parent population as the SDSS ellipticals. 

Our first paper (Haan et al. 2011) contains fitting results for 19 GOALS ellipticals 
presented in this paper. Haan et al. allow the Sersic index to vary. To be consistent with 
the present I-band analysis, we refit the 19 galaxies with GALFIT using the criteria outlined 
in Section 3. The result is presented in Table 8 and i\/2 vs. < fi\/2 > plot for these ellipticals 
in I-band (red circles) and in H-band (brown squares using I — H = 2.65) are presented in 
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Fig. 16b where the same object is connected with a dotted line. 

In general, the r 1( / 2 estimated from H-band is about 2.4 kpc smaller than that from 
I-band (medians for Ty 2 (I) and ry 2 (H) are 4.8 kpc and 2.4 kpc respectively), but no 
correlation is found between the ri/ 2 (I) and ri/ 2 (H) as shown in Fig. 17a. In contrast, a 
weak correlation was found between < /ix/ 2 > (/) and < //i/ 2 > (H) (Fig. 17b, r=0.46). 
We also find that the difference in half-light radius at I and H bands has a good correlation 
with the I-H color of the galaxy; the difference in the half-light radius become smaller for 
the objects with red I-H color (Fig. 17c, r=0.75). This suggests that the objects with blue 
I-H color tend to have strong nuclear emission at H that in turn causes a large difference 
between r^ (I) and r 1 / 2 (H). As shown in Fig. 17d, a strong correlation exists between 
the host magnitude of m T and m H (r=0.90). The mean (median) of mj - m H in the plot 
is 1.65±0.63 (1.78) mag. This value is slightly larger than that found in SDSS ellipticals 
(I — H = 1.41 mag, Jahnke et al. 2004) probably because the GOALS ellipticals have more 
dust and are thus redder. 



5.3. Comparison with Prior Work 

In an attempt to study the evolutionary connection between ULIRGs and QSOs (QUEST: 
Quasar ULIRG Evolution Study) Veilleux et al. (2006, 2009a) have performed morphological 
analysis of ULIRGs and PG QSO hosts using HST NICMOS H-band images. The QUEST 
ULIRG sample consists of 26 highly nucleated (i.e., the ratio of the H-band luminosity from 
the inner 4 kpc to the total luminosity is larger than 1/3) ULIRGs drawn from 1 Jy sample 
of ULIRGs (Kim & Sanders 1998). The majority of these galaxies host strong optical and 
infrared AGNs and are therefore AGN-biased and not a representative sample of the 1 Jy 
ULIRGs. 

One of the main distinctions of morphological types between the GOALS sample and 
QUEST sample is that elliptical fraction and spiral fraction in QUEST sample are larger than 
50% (74% and 54% in QUEST ULIRGs and PG QSO hosts, see Table 3) and less than 20% 
(16% and 4% in ULIRGs and PG QSO hosts) respectively. In contrast, the elliptical fraction 
in GOALS (U)LIRGs is only 5% when they are widely-separated and ~20% when they 
become single-nucleus mergers. Most of the QUEST galaxies are single-nucleus advanced 
mergers, and if all of the LIRGs evolve to ULIRGs, a significant fraction of the LIRGs turned 
into ellipticals only when they are completely merged. 

Interestingly, despite the large difference in morphological types between GOALS and 
QUEST samples, we do not see any difference in host galaxy magnitudes. The I-band 
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host absolute magnitude of the QUEST ULIRGs is -22.48±0.65 which is similar to that of 
the GOALS galaxies. A K-S test shows that the host absolute magnitudes of the GOALS 
galaxies are similar (P (null) =0.42) to that of the QUEST ULIRGs. The distribution of ab- 
solute magnitudes of the PG QSO hosts is also not much different from that of the GOALS 
galaxies. The host absolute magnitudes for the PG QSOs range from -21.54 to -24.43 mag 
with mean value of -22.87±0.78 mag. A K-S test of the absolute magnitudes between the 
GOALS galaxies and PG QSO hosts shows that both samples are not significantly differ- 
ent (P(null)=0.39). The negligible differences in the host galaxy magnitudes between the 
GOALS galaxies, ULIRGs, and PG QSOs suggests that the GOALS galaxies satisfy one of 
the necessary conditions to evolve into ULIRGs & QSO hosts. 

The PSF fraction in GOALS single-nucleus advanced mergers is 38%. The PSF frac- 
tion in QUEST ULIRGs is 100% probably because they were selected to be nucleated and 
therefore more likely to have dominant AGN. The mean (median) IpsF/^host of the QUEST 
ULIRGs and PG QSO hosts are 0.69±1.56 (0.10) and 3.52±4.58 (1.88) respectively. The 
median of the IpsF/lhost of the QUEST ULIRGs and PG QSO hosts are about 8 and 140 
times larger than that in GOALS galaxies. This can be caused by distance effects since the 
GOALS galaxies are much closer than QUEST ULIRGs and PG QSOs. The IpsF/^host as a 
function of redshift is plotted in Fig. 18 where red circles, blue circles, and blue filled circles 
represent GOALS galaxies, QUEST ULIRGs, and PG QSO hosts respectively. We find no 
distance dependency of Ipsp/^-host in GOALS and QUEST ULIRGs, except for a weak cor- 
relation observed in PG QSO hosts. The extinction in /-band is about 3 times larger than 
that in H-band. Thus, even if extinction is considered, the IpsF/^-host in GOALS galaxies 
and QUEST ULIRGs are significantly smaller than that in PG QSO hosts. 

Fig. 19 shows IpsF/^host as a function of /25//60 (Fig- 19a) and infrared luminosity (Fig. 
19b) where red circles, red circles with blue dot, blue circles, and blue filled circles represent 
GOALS LIRGs (I-band), GOALS ULIRGs (I-band), QUEST ULIRGs, and PG QSO hosts, 
respectively. The arrows represent either upper or lower limits and the circles with small, 
medium, and large crosses representing ULIRGs with outflow velocities less than 1000 km 
s- 1 (small crosses, F00188-0856, F03250+1606, F04313-1649, F05189-2524, F09039+0503, 
F09539+0857, F11506+1331, F14197+0813, F20414-1651; Rupke et al. 2005a), 1000 to 
5000 km s" 1 (medium cross, F05024-1941; Rupke et al. 2005b), and 5,000 to 10,000 km 
s- 1 (large crosses, F07599+6508, F12540+5708; Rupke et al. 2005b), respectively. If we 
consider QUEST ULIRGs and PG QSO hosts only, we find a strong correlation (solid line 
in Fig. 19a, r=0.78) between /25//60 and IpsF/^host- The IRAS color ratio of /25//60 is a 
well known indicator of AGN activity in infrared galaxies (e.g., de Grijp et al. 1985, Miley 
et al. 1985). Thus, this correlation could indicate that the PSF strength is directly related 
to the AGN activity. For a given /25//60, the IpsF/lhost in GOALS galaxies estimated from 
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I-band is about an order of magnitude smaller than that in QUEST ULIRGs. This can be 
explained by a combined effect of larger extinction at /-band (GOALS galaxies) than at 
i/-band (QUEST ULIRGs) and/or a relatively weak AGN strength in GOALS galaxies than 
in QUEST ULIRGs. 

In Fig. 19b, if only the GOALS galaxies and QUEST ULIRGs are considered, a weak 
correlation can be found between the and IpsF/^-host (r=0.48). However, if the PG QSO 
hosts are included, the correlation disappears. The PSF strength in QUEST ULIRGs tend 
to be proportional to the outflow velocity. The object with the largest outflow velocity 
has the strongest PSF strength. The intense starburst in a ULIRG will inevitably power 
strong superwinds which in turn will clear out gas and dust in the nuclear region. If such 
an event occurs, the spectral energy distribution in a ULIRG no longer peaks in the far- 
IR wavelength region; the peak instead shifts to optical/UV wavelengths. As a result, the 
ULIRGs will become proto-quasars (Sanders et al. 1988; Veilleux et al. 2006, 2009a, 2009b). 

The mean (median) l res /l host in GOALS galaxies, QUEST ULIRGs, and PG QSO hosts 
are 0.38±0.11 (0.36), 0.20±0.09 (0.20), and 0.14±0.06 (0.13) respectively. There is a clear 
tendency for decreasing l res /lhost from GOALS galaxies to QUEST ULIRGs and QUEST 
ULIRGs to PG QSO hosts. The l res /l host in GOALS ULIRGs is about twice that of QUEST 
ULIRGs. Compared to the //-band images, the /-band images are more sensitive to the 
detection of tidal features. In addition, the GOALS ULIRGs (mean NS=4.0±9.6 kpc) are 
less evolved than the QUEST ULIRGs (mean NS=1.0±4.0 kpc). The combined effect might 
explain the l re s/^-host difference between the GOALS ULIRGs and QUEST ULIRGs. The 
average l res /lhost in PG QSO hosts is about 30% less than that in QUEST ULIRGs. This 
indicates that the PG QSO hosts are in general more evolutionary advanced systems than 
the QUEST ULIRGs. 



6. Summary 

Hubble Space Telescope ACS/WFC F814W-band imaging data has been used to carry 
out a 2-dimensional structural analysis of 85 Lir > 10 1L4 L (U)LIRGs in GOALS sample. 
The main results of the study are as follows: 

1. The fraction of GOALS galaxies best fit by disk, disk+bulge, and elliptical profiles is 
32.1%, 48.9%, and 19.0%, respectively. The mean host absolute magnitude of the GOALS 
galaxies is — 22.64±0.62 mag, which corresponds to 1.8^q;4 L}. We do not find any difference 
of the mean host absolute magnitudes between the GOALS LIRGs and ULIRGs. 

2. The mean bulge absolute magnitude of GOALS galaxies is about 1.1 magnitude 
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fainter than the mean host magnitude. The median B/D for the GOALS disk+bulge sys- 
tems is 0.67 which corresponds to that of the Sa-Sab galaxies. Mass ratios in the GOALS 
binary systems suggest that most of the galaxies (93%) are the result of major mergers 

ij^ secondary / ^primary ^* 1/4)- 

3. The average bar fraction in GOALS galaxies (26%) is about one third of that found 
in nearby spirals (~65-70%). The bar fraction in widely separated (NS=50 kpc) GOALS 
galaxies is about 56%, but may be as high as 72% if orientation effects are taken into 
account. We find twice as many bars in GOALS disk+bulge systems (32.8%) than in pure- 
disk galaxies (15.9%). But most of the disk+bulge systems that contain bars are disk- 
dominated, suggesting that bars in disk+bulge systems are preferably found in large disks 
with small bulges. This suggests that large bulges in merging galaxies may be stabilizing 
the disk against bar formation. But note that this trend is different from the bar fraction in 
nearby normal, bright disk galaxies (bar fractions are about the same across Sa-Sc types). 
This trend is also opposite of the behavior seen at high redshifts where bars and bulges are 
coeval at z~l - bars are very scarce in very late type galaxies at high redshifts but evolve 
the fastest from z~l to the present to reach parity with early type systems. 

The physical properties of bars in GOALS galaxies such as the bar-to-host intensity 
ratio (1^/1^0^=12.8+7.8%), bar half-light radius (r e (bar)=2.7±1.6 kpc), and bar ellipticity 
(ebar=0.68±0.14) are similar to those found in nearby spiral galaxies. We do not find any 
trend of Ibar/^-host and r e (bar) as a function of nuclear separation, which may be an indication 
that bar destruction in merging galaxies occurs rather rapidly. 

4. A nuclear PSF component has been detected in 16% of the GOALS galaxies, with 
an average PSF-to-host intensity ratio of ^5^/1^4=0.01. The fraction of detected nuclear 
PSFs increases toward later merger stages and is highest in systems with a nuclear separation 
of kpc. Thus, light from an AGN and/or compact nuclear star cluster is more visible at 
I-band or AGN activity is more frequent as (U)LIRGs enter their latter stages of evolution. 

5. The GOALS galaxies have significant fraction of their light in tidal features (l re s/^-host 
= 0.38+0.11). Examination of the residual strength in GOALS binary systems suggests small 
companions suffer more tidal impact than larger companions (about 50% more). We do not 
observe a rapid decrease in tidal features in GOALS galaxies as a function of decreasing 
nuclear separation, suggesting a significantly long time is required to settling down the tidal 
features even when the progenitor galaxies are completely merged. The late-stage mergers 
whose profiles are well fit by a Sersic index of n=l may be in the process of reforming disks 
from infall. 



6. The mean size (^2=5. 3+3. 7 kpc) and slope of the surface brightness vs. half-light 
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radius relation (i.e., the Kormendy Relation) in GOALS elliptical hosts are similar to those 
in nearby SDSS ellipticals. The above results are consistent with the possibility that the 
GOALS galaxies belong to the same parent population as the SDSS ellipticals. 
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Table 1. Basic Properties of the GOALS Sample 
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11.71 


UGC 08696 


5 


0.0 


0.10 


173.0 


12.21 


MCG -03-04-014 





0.0 


0.12 


144.0 


11.65 


NGC 5331 


2 


17.4 


0.10 


155.0 


11.66 


CGCG 436-030 


2 


33.0 


0.14 


134.0 


11.69 


IRAS F14348-1447 


4 


6.2 


0.08 


387.0 


12.39 


IRAS F01364-1042 


5 


0.0 


0.07 


210.0 


11.85 


IRAS F14378-3651 


6 


0.0 


0.10 


315.0 


12.23 


III Zw 035 


3 


4.7 


0.08 


119.0 


11.64 


VV 340a 


1 


27.9 


0.06 


157.0 


11.74 


NGC 0695 





0.0 


0.11 


139.0 


11.68 


VV 705 


4 


6.9 


0.16 


183.0 


11.92 


MCG +05-06-036 


1 


32.4 


0.12 


145.0 


11.64 


ESO 099-G004 


3 


3.7 


0.14 


137.0 


11.74 


UGC 02369 


2 


13.5 


0.23 


136.0 


11.67 


IRAS F15250+3608 


5 


0.0 


0.18 


254.0 


12.08 


IRAS F03359+1523 


3 


7.7 


0.11 


152.0 


11.55 


UGC 09913 


5 


0.0 


0.08 


87.9 


12.28 


ESO 550-IG 025 


2 


11.4 


0.19 


138.5 


11.51 


NGC 6090 


4 


4.4 


0.19 


137.0 


11.58 


NGC 1614 


5 


0.0 


0.23 


67.8 


11.65 


IRAS F16164-0746 


5 


0.0 


0.06 


128.0 


11.62 


ESO 203-IG001 


4 


8.1 


0.08 


235.0 


11.86 


ESO 069-IG006 


2 


65.1 


0.10 


212.0 


11.98 


VII Zw 031 





0.0 


0.11 


240.0 


11.99 


IRAS F16399-0937 


3 


5.3 


0.13 


128.0 


11.63 


IRAS F05189-2524 


6 


0.0 


0.26 


187.0 


12.16 


NGC 6240 


4 


0.8 


0.15 


116.0 


11.93 


MCG +08-11-002 


6 


0.0 


0.08 


83.7 


11.46 


IRAS F17132+5313 


2 


6.6 


0.11 


232.0 


11.96 


IRAS F06076-2139 


3 


6.0 


0.10 


165.0 


11.65 


IRAS F17138-1017 





0.0 


0.14 


84.0 


11.49 


ESO 255-IG007 


2 


10.7 


0.16 


173.0 


11.90 


IRAS F17207-0014 


5 


0.0 


0.05 


198.0 


12.46 


AM 0702-601 


1 


56.7 


0.20 


141.0 


11.64 


IRAS 18090+0130 


2 


46.1 


0.10 


134.0 


11.65 


IRAS 07251-0248 


5 


0.0 


0.10 


400.0 


12.39 


IC 4687 


2 


31.7 


0.18 


81.9 


11.62 


IRAS 08355-4944 


3 


1.5 


0.24 


118.0 


11.62 


IRAS F18293-3413 


1 


5.6 


0.11 


86.0 


11.88 


NGC 2623 


5 


0.0 


0.08 


84.1 


11.60 


NGC 6670 


2 


17.0 


0.12 


129.5 


11.65 


ESO 060-IG 016 


3 


8.2 


0.13 


210.0 


11.82 


VV 414 


1 


38.0 


0.19 


113.0 


11.49 


IRAS F08572+3915 


4 


6.0 


0.24 


264.0 


12.16 


ESO 593-IG008 


2 


2.4 


0.08 


222.0 


11.93 


IRAS 09022-3615 


5 


0.0 


0.10 


271.0 


12.31 


IRAS F19297-0406 


5 


0.0 


0.09 


395.0 


12.45 


IRAS F09111-1007 


1 


39.2 


0.11 


246.0 


12.06 


IRAS 19542+1110 





0.0 


0.12 


295.0 


12.12 


UGC 04881 


2 


9.5 


0.10 


178.0 


11.74 


IRAS 20351+2521 





0.0 


0.12 


151.0 


11.61 


UGC 05101 


5 


0.0 


0.09 


177.0 


12.01 


CGCG 448-020 


2 


7.9 


0.18 


161.0 


11.94 


ESO 374-IG 032 


4 


2.7 


0.12 


148.5 


11.73 


ESO 286-IG019 


5 


0.0 


0.15 


193.0 


12.06 


IRAS F10173+0828 





0.0 


0.10 


224.0 


11.86 


IRAS 21101+5810 


2 


7.5 


0.10 


174.0 


11.81 


NGC 3256 


5 


0.0 


0.15 


38.9 


11.64 


ESO 239-IG002 


5 


0.0 


0.16 


191.0 


11.84 


IRAS F10565+2448 


2 


22.9 


0.10 


197.0 


12.08 


IRAS F22491-1808 


5 


0.0 


0.10 


351.0 


12.20 


MCG +07-23-019 


2 


11.0 


0.11 


158.0 


11.62 


NGC 7469 


2 


26.3 


0.22 


70.8 


11.65 


IRAS F11231+1456 


1 


50.8 


0.10 


157.0 


11.64 


ESO 148-IG002 


4 


4.2 


0.15 


199.0 


12.06 


NGC 3690 


3 


5.0 


0.22 


50.7 


11.93 


IC 5298 





0.0 


0.21 


119.0 


11.60 


IRAS F12112+0305 


4 


4.4 


0.08 


340.0 


12.36 


ESO 077-IG014 


2 


14.2 


0.08 


186.0 


11.76 


IRAS 12116-5615 





0.0 


0.12 


128.0 


11.65 


NGC 7674 


1 


18.9 


0.36 


125.0 


11.56 


CGCG 043-099 


5 


0.0 


0.09 


175.0 


11.68 


IRAS F23365+3604 


5 


0.0 


0.13 


287.0 


12.20 


ESO 507-G070 


6 


0.0 


0.06 


106.0 


11.56 


IRAS 23436+5257 


4 


3.3 


0.13 


149.0 


11.57 


IRAS 13120-5453 


5 


0.0 


0.07 


144.0 


12.32 


MRK 0331 


1 


43.6 


0.14 


79.3 


11.50 


VV 250a 


2 


22.1 


0.17 


142.0 


11.81 















Col 1: Object name. See Armus et al. (2009) for details 
Col 2: Interaction class (Haan et al. 2011, Evans et al. 2013) 



Col 3: Nuclear separation in kpc unit 
Col 4: IRAS /25 to /go flux ratio 
Col 5: Luminosity distance in mpc unit 
Col 6: Infrared luminosity of the object 
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Table 2. GALFIT Results 



IV (XL1H-, 
(1) 


(2) 


m n 
(3) 


r i 

2 

(4) 


Ml 

2 

(5) 


b 

a 

(6) 


m m 
(7) 


(8) 


J3 

(9) 


Iter 

thost 

(10) 


-W 

(11) 


i rcs 
(12) 


IVlf 

(13) 


Mi 
(14) 


NOT 00^4 


4 


12 


82 


2 


7 


18 


99 





87 


12 


42 


18.13 


N 






0.4 


35 


3 


-22.39 


-22.39 




4 


13 


69 





5 


16 


25 





69 
























A rn okc ivr 
r\.L \j ZOO IN 


4 


13 


70 


10 


3 


22 


11 





35 


13 


46 




Y 


15 


9 




84 





-22 83 


-22 83 




13 


15 


23 


2 


8 


20 


82 





37 
























Arn 2^6 S 

illJJ iitJVJ kJ 




13 


68 


3 


4 


19 


69 





61 


13 


56 




N 








56 











4 


16 


00 





7 


18 


41 





32 
























mco +12-02-001 s 


4 


14 


13 


3 


3 


21 


15 





32 


13 


55 


17.17 


Y 






2.9 


45 


7 


-21.31 


-21.29 




4 


14 


51 


2 


7 


21 


06 





48 
























unn 4-19 09 nni n 


4 


14 


24 


6 


2 


22 


61 





68 
















55 


5 






TO 1 623 w 


4 


13 


12 


1 


4 


17 


89 





92 


13 


12 




N 








42 


1 


-22.71 


-22.71 


TC 1 623 F, 


4 


12 


94 


5 





20 


45 





39 








N 








40 


9 






mcc -03-n4-ni4 

XV X V_ V.. I WJ u J- i 


4 


13 


52 


3 


1 


18 


89 





85 


13 


05 




N 








26 


8 


-22.79 


-22.79 




4 


14 


18 


3 


9 


20 


04 





57 
























mm 436 030 f, 


4 


15 


87 


7 


4 


23 


29 





51 


15 


51 




Y 


18 


7 




54 


5 


-22.20 


-22.19 




4 


18 


10 


2 


9 


23 


50 





31 


























13 


17 


30 


2 


6 


22 


47 





39 
























mm 436 030 w 

v. v. i v_ v.i i mnwu vv 


4 


13 


75 


4 


2 


19 


96 





64 


13 


75 


18.64 


N 






1.0 


45 


3 






TRAS F01 364-1 042 


4 


15 


17 


3 


5 


20 


03 





55 


15 


17 




N 








36 


6 


-21.54 


-21.54 


TTT Zw 035 N 


4 


15 


11 


1 


8 


19 


73 





57 


14 


04 




N 








18 





-21.47 


-21.47 




4 


14 


55 


1 


7 


19 


06 





38 
























TTT Zw 035 s 


4 


16 


36 


1 


1 


19 


90 





68 








N 








24 


4 






mpp 06Q5 


4 


12 


81 


4 


8 


19 


19 





95 


12 


39 




N 








38 


7 


-23.25 


-23.25 




4 


13 


63 


3 


4 


19 


29 





83 
























MCO 4-05-06-036 F 


4 


13 


34 


5 





19 


75 





62 


12 


93 




Y 


12 


5 




22 


5 


-23.48 


-23.47 




4 


14 


54 


1 


2 


17 


84 





50 


























13 


15 


58 


2 


8 


20 


70 





25 
























Men 4-0^ nfi n^fi w 

1 VI V_/ v_t yU J UU UOU VV 


4 


14 


56 


5 


5 


21 


17 





61 


13 


58 


18.84 


Y 


23 


3 


0.6 


25 


1 








4 


14 


49 


2 


1 


18 


97 





65 


























13 


15 


58 


2 


8 


20 


70 





25 
























UbU Uzooy IN 


i 
i 


13 


72 


7 


4 


21 


10 





34 


12 


85 




M 

IN 








18 


9 


-zo.Zz 


-zo.zl 




4 


13 


50 


4 


6 


19 


83 





78 
























UGC 02369 S 


1 


14 


12 


4 


5 


20 


41 





75 


13 


72 


17.58 


N 






2.6 


38 


7 








4 


14 


99 


3 


4 


20 


68 





71 
























IRAS F03359+1523 E 


1 


15 


08 


2 


9 


20 


18 





24 


15 


08 




N 








41 


3 


-21.53 


-21.53 


IRAS F03359+1523 W 


4 


15 


61 


1 


2 


18 


76 





87 
















5 


5 






ESO 550-IG 025 N 


1 


13 


91 


4 


6 


20 


21 





57 


13 


34 




N 








24 


9 


-22.77 


-22.77 




4 


14 


31 


3 


6 


20 


07 





67 
























ESO 550-IG 025 S 


1 


14 


90 


2 


8 


20 


13 





34 


14 


66 




N 








16 


8 








4 


16 


44 





9 


19 


29 





44 
























NGC 1614 


4 


11 


83 


5 


5 


20 


00 





80 


11 


83 


16.35 








1.5 


55 


5 


-22.36 


-22.34 


ESO 203-IG001 E 


1 


16 


72 


1 


6 


19 


62 





49 


15 


05 




N 








19 


2 


-21.89 


-21.89 




4 


15 


31 


9 





22 


01 





61 
























ESO 203-IG001 W 


4 


17 


05 


2 


5 


20 


96 





48 


17 


05 












35 









VII Zw 031 


1 


14 


25 


3 


4 


18 


78 





92 


14 


10 


19.65 


N 






0.5 


27 


3 


-22.99 


-22.99 




4 


16 


29 


1 





18 


22 





80 
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Table 2 — Continued 



(1) 


(2) 


m n 
(3) 


r i 

2 

(4) 


(5) 


6 

a 

(6) 


m m 
(7) 


m psf 

(8) 


13 

(9) 


har 
I host 

(10) 


Ipsf 
1 host 
(11) 


Ires 
Ihost 

(12) 


M, 
(13) 


aL host 

(14) 


TRAS F051 89-2524 


4 


14 


81 


8 


6 


21 


86 





85 


13 


66 


18.99 


N 




0.7 


35 


6 


-22.81 


-22.80 




4 


14 


12 


1 





16 


45 





93 
























1 


12 


92 


4 


1 


20 


03 





75 


12 


92 




N 






30 


2 


-21.90 


-21.90 


TRAS F06076-2139 N 

iiiriij i_ w w w 1 w iiuj i> 


I 


14 


28 


5 


1 


20 


47 





78 


13 


90 




Y 






16 


8 


-22.45 


-22.45 




4 


15 


22 


1 


9 


19 


25 





54 






















TRAS F06076-2139 S 


1 


15 


82 


6 


3 


22 


45 





70 


15 


19 




Y 


5.5 




22 


1 








4 


16 


28 


1 


1 


19 


04 





82 
























B 


18 


07 


3 


8 


23 


60 





12 






















FSO 255-TG007 NW 


4 


13 


73 


4 


6 


19 


58 





80 


13 


73 










34 


4 


-23.10 


-23.10 


fso 255-TO007 n 


4 


15 


22 


2 


1 


19 


35 





38 


14 


55 




Y 






26 


6 








4 


15 


39 


1 


8 


19 


18 





53 






















fso 2^5-10007 s 


4 


15 


09 


4 


4 


20 


84 





27 








N 






32 


5 






AM 0702-601 N 


4 


14 


21 


3 


4 


19 


84 





90 


14 


09 


18.22 


Y 




2.0 


33 


1 


-22.74 


-22.73 




4 


16 


50 





3 


16 


56 





64 






















AM 0702-601 S 


4 


13 


37 


4 


8 


19 


72 





87 














36 


6 






TRAS 07251-0248 


1 


17 


54 


1 


4 


19 


16 





65 


15 


40 




N 






41 


3 


-22.72 


-22.72 




4 


15 


56 


6 


4 


20 


52 





57 






















IRAS 08355-4944 E 


4 


14 


70 


4 


1 


21 


10 





62 


14 


70 


17.81 






5.6 


31 


3 


-21.03 


-20.98 


TRAS 08355-4944 W 


4 


15 


78 


2 





20 


59 





39 








N 






41 


3 






Nnn 2623 




13 


64 


2 


2 


19 


38 





57 


12 


33 




N 






38 


4 


-22.15 


-22.15 




4 


12 


71 


13 


6 


22 


40 





56 






















FSO 060-TO 01 fi E 


1 


15 


36 


3 


6 


20 


31 





34 


15 


36 




N 






25 


8 


-22.55 


-22.55 


FSO OfiO-TO 016 W 

1 V../ V ' KJ V 7 1 VI V / 1 \t vv 


4 


14 


60 


10 


3 


21 


82 





16 








N 






36 


6 






TRAS F085724-3915 SF 




16 


18 


6 


7 


22 


02 





36 


16 


08 




Y 


9.7 




67 


9 


-21.73 


-21.73 




13 


18 


73 





4 


18 


38 





50 






















TRAS F085724-3915 NW 

iiinij i_ woo i iii^ijt/iij i> vv 


4 


16 


39 


2 


6 


20 


14 





72 








N 






49 


2 






TRAS 00022-3615 


4 


14 


36 


6 


1 


19 


96 





82 


14 


36 










53 





-22.81 


-22.80 


TR AS F0Q1 1 1-1 007 F 


4 


14 


24 


7 


8 


20 


54 





66 


14 


04 




Y 


17.0 




29 


1 


-23.38 


-23.38 




B 


16 


68 


2 


7 


20 


72 





31 
























13 


16 


78 





2 


15 


49 





55 






















IRAS F09111-1007 W 


1 


15 


42 


4 


4 


20 


49 





70 


14 


72 




Y 






27 











4 


15 


53 


3 


2 


19 


91 





50 






















UGC 04881 NE 


1 


13 


65 


7 


1 


20 


40 





57 


13 


29 




Y 






44 


5 


-23.32 


-23.32 




4 


14 


68 


6 


9 


21 


37 





42 






















UGC 04881 SW 


1 


14 


96 


2 


5 


19 


40 





56 








N 






39 


8 






UGC 05101 


1 


14 


28 


4 


8 


20 


17 





87 


13 


47 


19.21 


N 




0.5 


28 


6 


-22.89 


-22.88 




4 


14 


16 


6 


2 


20 


60 





38 






















ESO 374-IG 032 W 


1 


15 


60 


1 


5 


19 


38 





74 


13 


22 




N 






42 


1 


-22.49 


-22.49 


ESO 374-IG 032 E 


4 


13 


35 


9 


3 


21 


08 





71 






















IRAS F10173+0828 


1 


16 


00 


4 


3 


21 


22 





17 


15 


13 




N 






26 


3 


-21.64 


-21.64 




4 


15 


77 


3 


5 


20 


52 





35 






















NGC 3256 


1 


10 


79 


3 


1 


18 


89 





69 


10 


70 




N 






47 


4 


-22.59 


-22.59 




4 


13 


46 





3 


16 


81 





67 






















IRAS F10565+2448 W 


4 


13 


76 


4 


8 


19 


47 





78 


13 


76 










33 


4 


-23.00 


-23.00 


IRAS F10565+2448 E 


4 


15 


38 


1 


3 


18 


24 





90 














14 


6 
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Table 2 — Continued 



Name 


n m n 


r i 

2 




6 
a 




B 


^bar 
I host 


-w 


J re5 


Mt M host 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 


(11) 


(12) 


(13) (14) 



IRAS F10565+2448 S 


4 


17 


13 





9 


19 


08 





74 












13 


3 






MCG +07-23-019 


1 


13 


94 


4 


4 


19 


90 





25 


13.94 




N 






56 


5 


-22.43 


-22.43 


IRAS F11231+1456 W 


1 


13 


67 


9 





21 


17 





34 


13.31 




Y 


17.4 




37 


3 


-23.13 


-23.13 




4 


15 


95 


1 





18 


61 





65 






















B 


15 


11 


6 


1 


21 


78 





12 




















IRAS F11231+1456 E 


1 


14 


44 


3 


9 


20 


11 





29 






Y 






37 


3 






NGC 3690 SW 


1 


11 


86 


6 


9 


21 


14 





48 


11.32 


16.99 


N 




0.7 


45 


3 


-22.79 


-22.78 




4 


12 


34 


2 


4 


19 


37 





55 




















NGC 3690 NE 


1 


11 


95 


2 


6 


19 


10 





83 






N 






33 


4 






IRAS F12112+0305 N 


1 


16 


11 


3 


3 


19 


91 





64 


16.11 




N 






49 


7 


-22.68 


-22.68 


IRAS F12112+0305 S 


4 


15 


64 


9 


3 


21 


70 





84 




















IRAS 12116-5615 


4 


13 


77 


4 





19 


93 





94 


13.77 


18.02 






2.0 


22 


1 


-21.78 


-21.75 


CGCG 043-099 


1 


15 


33 


2 


9 


20 


17 





61 


13.45 




Y 






33 


4 


-22.64 


-22.64 




4 


13 


66 


11 


9 


21 


56 





74 




















ESO 507-G070 


1 


13 


36 


8 


2 


21 


47 





30 


12.36 




N 






28 


6 


-22.58 


-22.58 




4 


12 


91 


18 


8 


22 


82 





41 




















IRAS 13120-5453 


1 


13 


86 


7 


9 


21 


27 





58 


13.38 


18.13 


N 




1.2 


32 


5 


-22.47 


-22.46 




4 


14 


51 


1 


2 


17 


89 





87 




















VV 250a E 


1 


16 


30 


1 


8 


20 


53 





54 


14.01 




Y 






46 


6 


-22.70 


-22.70 




4 


14 


15 


3 





19 


45 





51 




















VV 250a W 


1 


14 


71 


4 


1 


20 


70 





35 


14.14 


19.63 


N 




0.5 


50 


1 








4 


15 


11 


5 


6 


21 


80 





68 




















UGC 08387 


1 


13 


50 


4 


4 


20 


17 





47 


12.95 




N 






39 


8 


-22.48 


-22.48 




4 


13 


94 


10 


1 


22 


43 





62 




















NGC 5256 NE 


1 


13 


80 


4 


1 


20 


02 





69 


13.29 




N 






19 


8 


-23.37 


-23.37 




4 


14 


35 


1 


5 


18 


36 





59 




















NGC 5256 SW 


4 


12 


71 


8 


3 


20 


45 





79 












30 


8 






Arp 240 E 


1 


12 


41 


9 


2 


20 


73 





36 


12.04 




N 






53 


5 


-23.94 


-23.94 




4 


13 


39 


6 


5 


20 


96 





60 




















Arp 240 W 


1 


12 


16 


6 


8 


19 


81 





70 


12.14 




Y 


0.4 




55 





-23.94 


-23.94 




4 


17 


15 





2 


16 


87 





94 






















B 


17 


96 





4 


19 


20 





28 




















UGC 08696 


1 


14 


04 


3 


6 


19 


34 





54 


13.22 




N 






44 


1 


-23.17 


-23.17 




4 


13 


90 


10 





21 


44 





35 




















NGC 5331 N 


1 


13 


43 


5 


5 


19 


89 





50 


13.23 




Y 


4.9 




30 


5 


-23.66 


-23.66 




4 


15 


58 





4 


16 


47 





60 






















B 


16 


39 


2 





20 


69 





24 




















NGC 5331 S 


1 


13 


10 


7 


5 


20 


2d 





35 






N 






38 


4 






IRAS F14348-1447 NE 


1 


16 


07 


5 


3 


20 


67 





41 


16.07 




N 






29 


1 


-23.39 


-23.39 


IRAS F14348-1447 SW 


4 


14 


86 


14 





21 


57 





78 




















IRAS F14378-3651 


1 


16 


16 


4 


2 


20 


63 





67 


15.51 


19.88 


N 




1.7 


27 


5 


-22.03 


-22.01 




4 


16 


37 


1 


7 


18 


88 





90 




















VV 340a S 


1 


14 


01 


8 


2 


21 


31 





80 


13.24 




N 






26 


3 


-23.56 


-23.56 




4 


13 


98 


7 


5 


21 


10 





82 




















VV 340a N 


1 


13 


16 


8 


1 


20 


43 





37 






N 






35 


3 
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Table 2 — Continued 



1 N cXll 

(1) 


(2) 


m„ 
(3) 


r i 

2 

(4) 


(5) 


b 
a 

(6) 


m m 
(7) 


m psf 

(8) 


73 

(9) 


har 

(10) 


hsf 
Ihost 

(11) 


Ire.* 
I host 

(12) 


(13) 


wi -host 

(14) 


VV 705 N 


I 


14 


61 


3 


2 


19 


56 





78 


14 


33 






Y 






40 


9 


-22 


94 


-22.94 




4 


15 


95 





7 


17 


62 





67 


























VV 705 S 


[ 


14 


67 


4 


5 


20 


35 





43 










N 






57 


5 








F,SO 099-0004 N 


I 


14 


08 


8 


6 


21 


78 





51 


14 


08 






j\ 






38 





-22 


25 


-22.25 


f,so 099-0004 s 


I 


14 


42 


5 


3 


21 


06 





35 










N 
















TRAS FT52504-3608 


1 


15 


70 


1 


9 


18 


91 





69 


14 


88 


19 


38 


N 




1.3 


40 


5 


-22 


34 


-22.33 




4 


15 


57 


4 


5 


20 


63 





81 


























TJOO 0991 3 


1 


12 


85 


4 


3 


19 


93 





87 


11 


79 






N 






33 


7 


-22 


75 


-22.75 




4 


12 


31 


15 


5 


22 


20 





84 


























Mnn finqn w 

I'VJV VJVJ CVJ VV 


4 


15 


08 


1 


3 


18 


64 





42 


15 


08 






N 






81 


7 


-22 


81 


-22.81 


ngc 6090 e 


4 


13 


05 


4 


8 


19 


48 





88 
















43 


2 








TRAS FT 6T 64-0746 


I 


13 


79 


4 


2 


20 


05 





55 


13 


79 






N 






41 


7 


-22 


15 


-22.15 


FSO 069-TG006 N 


1 


14 


07 


8 


4 


20 


83 





46 


13 


81 






Y 


17.7 




53 


5 


-23 


65 


-23.65 




13 


15 


49 


2 


5 


19 


64 





42 


























fso 069-tgoo6 s 


1 


15 


40 


10 


9 


22 


73 





35 


14 


31 






N 






59 


2 










4 


14 


80 


2 


7 


19 


13 





82 


























TRAS FT 6399-0937 N 




14 


23 


4 


1 


20 


47 





70 


14 


23 






N 






45 


7 


-22 


16 


-22.16 


TRAS FT6399-0937 s 


4 


13 


67 


8 


2 


21 


41 





80 
















32 


8 








Non R240 




11 


75 


11 


1 


20 


34 





45 


11 


67 


18 


22 


N 




0.2 


39 


8 


-23 


79 


-23.79 




4 


14 


56 





9 


17 


67 





53 


























TRAS FT 7T 324-53T 3 F 




14 


97 


3 


8 


19 


84 





41 


14 


77 






Y 


13.7 




22 


5 


-22 


75 


-22.75 




B 


16 


70 


1 


6 


19 


69 





29 


























TRAS FT7T32+53T3 W 


1 


16 


63 


3 


7 


21 


42 





61 


15 


11 






N 






31 


6 










4 


15 


42 


2 


3 


19 


23 





83 


























TR AS FT 7T 38-1 (11 7 

a i in o ± j_ ( loo lui i 




13 


57 


3 


8 


20 


51 





59 


13 


24 






N 






23 


1 


-21 


27 


-21.27 




4 


14 


68 


5 


6 


22 


43 





46 


























IRAS F17207-0014 


1 


14 


35 


3 


7 


19 


48 





80 


14 


35 






N 






38 





-22 


43 


-22.43 


TRAS T8090+0130 F 




13 


93 


5 


4 


20 


66 





48 


13 


74 






N 






47 


9 


-22 


51 


-22.51 




4 


15 


75 


2 





20 


31 





38 


























td aq 1 enQn_L.ni ^n "W 


i 


14 


57 


5 


4 


21 


30 





51 


14 


49 






v 
i 


^ Q 




80 


9 










4 


18 


41 





5 


20 


05 





31 




























B 


17 


85 


2 


8 


23 


18 





17 


























IC 4687 N 


1 


12 


57 


3 





19 


05 





64 


12 


57 






N 






42 


9 


-22 


95 


-22.95 


IC 4687 W 


4 


13 


92 


1 


9 


19 


41 





67 
















34 


4 








IC 4687 S 


1 


13 


06 


3 


6 


19 


91 





52 


12 


88 






N 






25 


4 










4 


14 


92 


1 


3 


19 


57 





36 


























IRAS F18293-3413 NW 


4 


14 


56 





7 


17 


79 





61 


14 


56 












14 


3 


-22 


33 


-22.33 


IRAS F18293-3413 SE 


1 


12 


71 


2 


4 


18 


56 





64 










N 






38 











NGC 6670 E 


1 


16 


25 


2 





20 


93 





59 


16 


25 






N 






25 


1 


-22 


81 


-22.81 


NGC 6670 N 


1 


14 


60 


5 


2 


21 


30 





55 


13 


50 






N 






23 


3 










4 


13 


99 


6 


9 


21 


33 





20 


























NGC 6670 W 


1 


13 


80 


6 


5 


20 


99 





18 










N 






47 





-22 


81 


-22.81 


VV 414 NE 


1 


13 


82 


5 


6 


20 


96 





69 


13 


05 


17 


41 


N 




1.7 


54 





-23 


37 


-23.36 




4 


13 


78 


2 


5 


19 


21 





82 


























VV 414 SW 


1 


12 


99 


6 


6 


20 


51 





91 


12 


51 






Y 


17.0 




39 


8 
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Table 2 — Continued 



Name 


n 


m„ 


r I 


Mi 


b 

a 


m m 


m psf 


B 


'bar 
Ihost 


T 

Ihost 


I res 
Ihost 


M t 


M host 


(1) 


(2) 


(3) 


(4) 


(5 2 ) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 




4 


14.42 


1.6 


18.86 


0.84 




















13 


14.34 


3.7 


20.61 


0.26 


















ESO 593-IG008 W 


1 


14.04 


10.8 


21.26 


0.33 


14.04 




N 






20.3 


-23.19 


-23.19 


ESO 593-IG008 E 


1 


14.79 


7.1 


21.09 


0.20 






N 






98.2 






IRAS F19297-0406 


4 


15.38 


10.4 


21.42 


0.69 


15.38 


19.67 






1.8 


41.3 


-22.66 


-22.64 


IRAS 19542+1110 


4 


15.00 


1.3 


17.09 


0.92 


15.00 


19.59 






1.5 


21.7 


-22.34 


-22.32 


IRAS 20351+2521 


1 


13.14 


7.6 


20.36 


0.94 


12.85 




Y 


23.3 




34.4 


-23.04 


-23.04 




B 


14.77 


5.1 


21.11 


0.31 




















B 


15.89 


0.6 


17.62 


0.61 


















CGCG 448-020 W 


1 


13.86 


5.0 


20.06 


0.59 


13.86 




N 






48.8 


-22.79 


-22.79 


CGCG 448-020 S 


1 


14.82 


1.2 


17.84 


1.00 






N 






23.3 






CGCG 448-020 E 


1 


16.16 


3.3 


21.45 


0.79 






N 






25.8 






ESO 286-IG019 


4 


13.34 


7.8 


20.13 


0.78 


13.34 










61.4 


-22.97 


-22.97 


IRAS 21101+5810 NW 


4 


14.57 


9.1 


21.89 


0.67 


14.57 










43.6 


-21.65 


-21.65 


IRAS 21101+5810 SE 


4 


16.59 


1.0 


19.06 


0.76 












7.7 






ESO 239-IG002 


4 


12.71 


12.2 


20.49 


0.83 


12.71 


17.30 






1.7 


34.0 


-23.55 


-23.53 


IRAS F22491-1808 


1 


16.26 


4.5 


20.68 


0.28 


15.20 




N 






62.5 


-22.78 


-22.78 




4 


15.72 


4.9 


20.34 


0.44 


















NGC 7469 S 


1 


12.46 


5.6 


20.60 


0.53 


11.59 




Y 






36.3 


-23.13 


-23.13 




4 


12.24 


1.1 


16.88 


0.80 


















NGC 7469 N 


1 


12.89 


4.8 


20.70 


0.38 






N 






54.5 






ESO 148-IG002 N 


1 


14.65 


2.5 


18.90 


0.65 


13.94 




N 






49.2 


-22.97 


-22.97 


ESO 148-IG002 S 


4 


14.73 


3.8 


19.90 


0.75 


















IC 5298 


1 


14.33 


3.6 


20.43 


0.34 


12.79 




Y 


3.8 




18.5 


-22.60 


-22.55 




4 


13.15 


4.6 


19.77 


0.88 




















B 


16.39 


0.3 


17.09 


0.51 


















ESO 077-IG014 N 


1 


14.51 


4.8 


20.31 


0.36 


14.07 




Y 


15.2 




32.2 


-23.02 


-23.02 




4 


16.16 


2.0 


20.05 


0.84 




















B 


15.89 


1.4 


18.98 


0.64 


















ESO 077-IG014 S 


1 


15.01 


3.9 


20.38 


0.46 


14.59 




N 






42.9 








A 
1 


10.00 


A A 


01 A A 
Z1.44 


n ao 


















NGC 7674 W 


1 


12.43 


8.1 


20.17 


0.92 


12.23 




Y 


7.0 




30.5 


-23.50 


-23.50 




4 


14.70 


0.4 


16.10 


0.75 




















B 


15.12 


4.2 


21.43 


0.25 


















NGC 7674 E 


1 


14.87 


2.5 


20.06 


0.59 


13.89 




N 






31.6 








4 


14.45 


1.1 


17.94 


0.65 


















IRAS F23365+3604 


1 


14.59 


5.7 


19.92 


0.70 


14.51 




N 






38.4 


-22.99 


-22.99 




4 


17.33 


0.9 


18.53 


0.62 


















IRAS 23436+5257 N 


1 


14.24 


3.7 


19.91 


0.47 


14.07 




N 






52.0 


-22.34 


-22.34 




4 


16.19 


0.4 


16.96 


0.89 


















IRAS 23436+5257 S 


4 


15.19 


1.9 


19.47 


0.67 












18.9 






MRK 0331 E 


1 


14.44 


3.4 


21.25 


0.36 


12.43 




N 






29.1 


-22.25 


-22.25 




4 


12.62 


1.8 


18.03 


0.84 


















MRK 0331 W 


1 


14.49 


4.5 


21.91 


0.26 






N 






36.0 







Col 1: Object name 

Col 2: Best fit model components in each nucleus (1: Scrsic index n=l (disk), 4: Sersic index of n=4 (bulge), B: bar component) 
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Col 3: Apparent I-band magnitude for each component 
Col 4-' Half-light radius in kpc for each component 

Col 5: Mean surface brightness within half-light radius in I mag arcsec - 1 
Col 6: Axis ratio 

Col 7: Apparent I-band magnitude for all model components 
Col 8: Apparent I-band magnitude for PSF component 

Col 9: Visually-identified bar (Y). If no bar is visually identified, then a 'N' is placed in this column. If GALFIT can decompose 
bar, it is indicated as B in column (2) 

Col 10: Bar to host intensity ratio (%) 

Col 11: PSF to host intensity ratio (%) 

Col 12: Residual (total - model) to host intensity ratio (%) 

Col 13: Total absolute magnitude of the system 

Col 14-' Host absolute magnitude (total - PSF) of the system 
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Table 3. Distribution of Morphological Type 



Morphological Type 


GOALS 


GOALS LIRG 


GOALS ULIRG 


QUEST ULIRG 1 


PG QSO host 


Disk (%) 


32.1 


34.5 


22.2 


15.8 


3.8 


Disk+Bulge (%) 


48.9 


50.0 


44.4 


10.5 


42.3 


Elliptical (%) 


19.0 


15.5 


33.3 


73.7 


53.8 



1 : These ULIRGs were selected based on having a high degree of nucleation at H-band 



Table 4. Bar fraction 



Category 




Bar Fraction (%) 




N.S. (kpc) 





0-10 10-30 


30 - 70 




5.0 


17.9 41.7 


57.1 


Int. class 


1 


2 3 4 5 


6 




28.6 56.3 


30.8 17.7 16.7 6.3 


0.0 


Mor. type 


Disk 


Disk+Bulge 


Total 




15.9 


32.8 


25.9 



Table 5. Statistics of Barred and Unbarred Systems 



quantity 


M(host) mag 


M(bulge) mag 


B/D (median) 


D(%) 


D+B(%) 


barred 


-22.26±0.77 


-20.45±1.43 


0.36 


15.9 


32.8 


unbarred 


-22.03±0.80 


-21.23±0.88 


0.87 


84.1 


67.2 
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Table 6. PSF Detection Fraction 



Category 


PSF Detection Fraction (%) 




N.S. (kpc) 





0-10 10-30 


30 - 70 




38.5 


7.7 7.1 


17.4 


Int. class 


1 


2 3 4 5 


6 




33.3 16.7 


6.5 13.6 5.9 33.3 


50.0 


Mor. type 


Disk 


Disk+Bulge Elliptical 


Total 




2.3 


22.4 23.1 


16.1 


Infrared lum. 


LIRG 


ULIRG 






13.8 


25.0 





Table 7. Statistics of PSF and Non-PSF Systems 



quantity 


M(host) mag 


M (bulge) mag 


B/D 




PSF 


-22.30±0.71 


-21.34±1.00 


0.73±0.48 


0.16±0.06 


non-PSF 


-21.95±0.89 


-21.08±1.26 


1.37±1.83 


0.13±0.05 
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Table 8. I and H Bands Ty 2 and < //i/ 2 > for GOALS Galaxies with Eliptical Profiles 











I 












H 






Name 




IRpsF 


ri/2 


< Ml/2 > 




mpsF 


ri/2 


< Ml/2 > 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


IRAS F03359+1523 


15 


.61 




1.2 


18.76 


13 


.65 


17. 


.77 


1 


.4 


22.79 


NGC 1614 


11 


.83 


16.35 


5.5 


20.00 


10 


.90 


13 


.44 


0. 


.2 


17.94 


ESO 203-IG001 


1 7 
1 1 


.uo 




2.5 


20.96 




.oO 






q 
O 


A 

,4 


25.51 


ESO 255-IG007 


13 


.73 




4.6 


19.58 


11 


.62 






2 


,8 


21.97 


AM 0702-601 


13 


.37 




4.8 


19.72 


11 


.21 


18. 


.68 


4 


.1 


22.82 


IRAS F10565+2448 


13 


.76 




4.8 


19.47 


12 


.31 


16. 


.50 


1 


.5 


21.10 


IRAS F12112+0305 


15 


.64 




9.3 


21.70 


15 


.16 


16. 


.89 


1 


.3 


22.55 


IRAS 12116-5615 


13 


.77 


18.02 


4.0 


19.93 


10 


.93 


15. 


.46 


2 


.7 


21.82 


NGC 5256 


12 


.71 




8.3 


20.45 


11 


.48 






1 


.6 


23.65 


IRAS F14348-1447 


14 


.86 




14.0 


21.57 


14 


.17 


18. 


.99 


2. 


.1 


22.39 


NGC 6090 


13 


.05 




4.8 


19.48 


11 


.39 


19. 


.00 


5. 


.4 


23.63 


IRAS F16399-0937 


13 


.67 




8.2 


21.41 


11 


.89 


18. 


.08 


2 


.8 


22.88 


IRAS F19297-0406 


15 


.38 


19.67 


10.4 


21.42 


13 


.31 


17. 


.03 


5 


.5 


23.56 


IRAS 19542+1110 


15 


.00 


19.59 


1.3 


17.09 


13 


.02 


16. 


.03 


0. 


.6 


19.15 


ESO 286-IG019 


13 


.34 




7.8 


20.13 


12 


.08 






2 


.8 


22.24 
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Fig. 1. — Results from the GALFIT analysis. The upper panel shows the original image 
(a), model galaxy (b), and residual. The lower panel shows Sersic components used in the 
model. The component names listed on top of each plot are PSF, disk (n=l), and bulge 
(n=4). The intensity scale is logarithmic and stretch values are the same in all images. The 
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Fig. 2. — a) Distribution of morphological types in GOALS galaxies, GOALS LIRGs, and 
GOALS ULIRGs. The red, green, and blue colors represent disk, disk+bulge, and elliptical 
host respectively. 
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Fig. 3.- 



Distribution of morphological types as a function of interaction class. 
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Fig. 4. — Distribution of a) host absolute magnitude and b) bulge absolute magnitude for 
GOALS galaxies (red), GOALS LIRGs (green), and GOALS ULIRGs (blue). The vertical 
dashed line represents M\ = —22.0 mag, i.e., the /-band absolute magnitude of an L* galaxy 
in a Schechter luminosity function of the local field galaxies. 
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Fig. 5. — Distribution of bulge-to-disk ratio, B/D, in GOALS galaxies. 
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Fig. 6. — a) Distribution of mass ratio in binary systems and b) mass ratio as a function of 
M(host). 
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Fig. 7. — Bar fractions as a function of a) interaction class, and b) morphological type. 



-40- 




Fig. 8. — a) Distribution of Ibarfihost, and b) bar effective radius r e (bar). 
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Fig. 9. — a) Distribution of bar ellipticity, bar ellipticity as a function of b) bar half-light 
radius, and c) bulge-to-disk ratio. 
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Fig. 10. — PSF detection fractions as a function of a) interaction class, b) morphological 
type, and c) infrared luminosity (LIRG vs. ULIRG). 
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Fig. 11. — Distribution of IpsF/Ihost in GOALS galaxies (red), GOALS LIRGs (green), 
and GOALS ULIRGs (blue). The IpsF/^-host derived from I-band images is plotted on the 
left panel, and the ratio derived from H-band images for the GOALS galaxies classified as 
ellipticals is plotted on the right panel. The cyan line on top-left panel represents IpsF/^-host 
in GOALS elliptical galaxies. 
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Fig. 12. — a) IpsF/Ihost as a function of morphological type, b) M(bulge) vs. M(PSF) plot. 
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Fig. 13. — Distribution of l res /Ihost in GOALS galaxies (red), GOALS LIRGs (green), and 
GOALS ULIRGs (blue). 
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Fig. 14. — Host absolute magnitude vs. l res /lh os t for binary systems. In a binary system, if 
the lres/lhost i n a small companion (faint component) is larger than that in large one (bright 
component), it is plotted on the left panel (a) and the opposite case is plotted on the right 
panel (b). In both plots, the interacting pair is connected with a line, and green filled circles 
and red bars represent galaxies with PSF and stellar bar components, respectively. 
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Fig. 15. — a) Fraction of morphological types, b) host magnitude, c) huige/^-host, d) bar 
fraction, e) har/host, f) bar size, g) PSF fraction, h) IpsF/host, and i) l res /l host as a function 
of nuclear separation. 
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Fig. 16. — Surface brightness vs. half-light radius in the I-band for elliptical hosts in GOALS 
galaxies (red circles). The green and red lines represent least-square fits of the SDSS ellip- 
ticals and GOALS galaxies, respectively. The dotted lines in Fig. 16b connect the same 
objects observed in I-band (circles) and H-band (squares). 
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Fig. 17. — a) Half-light radius and b) surface brightness for GOALS elliptical hosts in I and 
H bands, c) I-H color vs. ri/ 2 (I)— ri/ 2 (H), and d) host magnitude of ni/ vs. m# plots. 
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Fig. 18. — PSF to host intensity ratio IpsF/^-host as a function of redshift for a) GOALS 
galaxies, b) QUEST ULIRGs, and c) PG QSOs. 
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Fig. 19. — IpsF/Ihost as a function a) /25//60, and b) Lir. The red circles, red circles with 
blue dot, blue circles, and blue filled circles represent GOALS LIRGs (I-band), GOALS 
ULIRGs (I-band), QUEST ULIRGs (H-band), and PG QSO hosts (H-band), respectively. 
The arrows represent either upper or lower limits and circles with small, medium, and large 
crosses represent QUEST ULIRGs with outflow velocity less than 1,000 km s _1 , 1,000 to 
5,000 km s _1 , and 5,000 to 10,000 km s _1 respectively. 
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